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Using lattice QCD, a diquark can be studied in a gauge-invariant manner by binding it to a static 
quark in a heavy-light-light hadron. We compute the simultaneous two-quark density of a diquark, 
including corrections for periodic boundary conditions. We define a correlation function to isolate 
the intrinsic correlations of the diquark and reduce the effects caused by the presence of the static 
quark. Away from the immediate vicinity of the static quark, the diquark has a consistent shape, 
with much stronger correlations seen in the good (scalar) diquark than in the bad (vector) diquark. 
We present results for nij^ = 293 MeV mNf — 2-\- \ QCD as well as nij^ — 940 MeV in quenched 
QCD, and discuss the dependence of the spatial size on the pion mass. 
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1. Introduction 

Diquarks are two-quark systems. Collective behavior of a diquark has been invoked to ex- 
plain many phenomena of strong interactions [1]. By introducing diquarks as effective degrees of 
freedom in chiral perturbation theory, they have been used to explain the enhancement of A/ = ^ 
nonleptonic weak decays [2]. A simple quark-diquark model is quite successful at organizing the 
spectrum of excited light baryon states [3] . 

The simplest diquark operators are quark bilinears with spinor part q^CTq. The favored com- 
binations are color antitriplet, even parity [4]. These are divided into "good" and "bad" diquarks. 
The good diquarks, q^Cy^q, have spin and are flavor antisymmetric due to fermion statistics. The 
bad diquarks, q^Cyiq, have spin 1 and are flavor symmetric. 

Both one-gluon exchange in a quark model [5, 6] and instanton [7] models give a spin coupling 
energy proportional to 5, • Sj, which favors the good diquark over the bad diquark. The strength 
of this coupling falls off with increasing quark masses. For the instanton model, the effective 
interaction has a flavor dependence that also favors the good diquark. 

2. Earlier studies in baryons 

Since diquarks are not color singlets, studying them within the framework of lattice QCD 
typically requires that they be combined with a third quark to form a color singlet. Diquark attrac- 
tions result in spatial congelations between the two quarks in the diquai'k, which can be probed by 
computing a wavefunction or two-quark density. 

In one study [8], using gauge fixing and 
three quarks of equal mass, a wavefunction was 
computed by displacing quarks at the sink, 

i/A(ri,r2)°<:^(M(r, + ri,?)c?(r, + r2,0^(r.,0 
xm(0,0)J"(0, 0)5(0,0)). 

Then, using the more convenient coordinates ^. , ^ ^ , , , , • , • 

ligure 1: Geometry or the three quarks used m this 

'fi=(r^+r2)/2 and r = (ri - r2)/2 (2 1) paper and in [8] (left), and restricted geometry used 

in [9] (right). 

(Fig. 1, left), the wavefunction of the good and 

bad diquarks was shown for different fixed /? = |R| as a function of r, in both Coulomb gauge and 
Landau gauge. In all cases, the wave function had a peak near r = 0, but it was found to fall off 
more rapidly for the good diquark, consistent with the expectation that good diquarks are more 
tightly bound. 

In a second study [9], spatial correlations were investigated by computing the two quark 
density P2{i^u,i^d) for a {u,d) diquark in the background of a static quai^k. To isolate correla- 
tions caused by the diquark interaction, analysis was restricted to spherical shells |ri,| = \rd\ = r 
(Fig. 1, right). For both good and bad diquarks, the density was found to be concentrated neai^ 
rud = |r„ — r^l =0, and the effect was much stronger for the good diquark. Fitting p2 for the 
good diquark to exp(— r„t//ro(r)), ro reached a plateau for lai^ge r, giving a characteristic size 
ro = l.l±0.2fm. 
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3. Correlation function 



In the first study, the wavefunctions of good and bad diquarks were compared for an unre- 
stricted geometry, but they were not compared against an uncorrelated wavefunction. In the sec- 
ond study, the intrinsic clustering caused by the diquark interaction was shown, however this was 
achieved by using a restricted geometry. 

To overcome these limitations, we combined a diquark with a static quark, using the baryon 
operator B = e"'"^ (mJCFJ^) s^, taking F = 75 for the good diquark and F = 71 for the bad diquark, 
and calculated the single quark density and the simultaneous two-quark density: 

n M = .V (0|g(0^^/yg('-,0g(0,O|0) . {0\B{0,tf)J-{r,,t)Jf,{r2,t)B{0jim 
P'^""^ {0\B{0,tf)B{(i,tm ' P'^'''''''' {0\BiO,tf)BiO,um 

Here, there are insertions of the current 7^ = /7^/, and the normalization factors A^i^2 are required 
since this local current is not conserved on the lattice. 

In a system where pi (r) is uniform, the two-particle correlation can be defined as 

Co(ri,r2) =p2(ri,r2) -pi(ri)pi(r2). 

Deviations from zero are seen as evidence for interactions between particles. This correlation 
integrates to zero and approaches zero as the relative distance ru = |ri — r2| increases beyond the 
range of interactions in the system. 

The situation considered here is not so simple. The single particle density is not uniform: it 
is concentrated near the static quark. Co will still integrate to zero and fall off at large distances, 
however it is also larger near the static quark and this obscures the diquark correlations. 

In order to remove the effect of the static quark, we define the normalized correlation function: 

r(r- r ^ P2(i-i.r2)-pi(ri)pi(r2) 

C(ri,r2j = . (3.1) 

Pi(ri)pi(r2) 

This divides out the tendency to stay near the static quark and retains the property of being zero if 
the two light quarks are uncorrelated (i.e. if p2(ri,r2) = pi(ri)pi(r2)). The downsides are that C 
no longer integrates to zero, and it is possible for C(r,r) to increase without bound as |r| — )• 00. 



4. Density in a periodic box 

We assume the lattice spacing is small enough that in an infinite volume we can treat p2 (ri , r2) 
as a function of/? = |R|,r=|r|, and 6 , the angle between R and r . Since the calculation is actually 
carried out on a finite lattice volume, to recover the infinite volume result, we need to deal with the 
effect of periodic boundary conditions. 

The problem of dealing with p in periodic boundary conditions has been previously analyzed 
for the case of a meson [10]. It was found that Pi(r) = £„g23 Pi (r + n^), where pi differs from the 
infinite volume result only for r > L due to interactions with periodic images. 

For this study, we have a baryon, and there is an additional complication: the contribution 
from "exchange diagrams" in which the two quarks travel in opposite directions across the periodic 
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boundary and can form a color singlet. As the lattice size grows, this becomes dominated by the 
propagation of the lightest meson and so falls off as exp(— m^rL). 

Ignoring the exchange diagrams and interactions with periodic images we find 

P2(ri,r2)= ^ P2(ri+iiiL,r2 + n2L), 

where P2 is the infinite volume two quark density. In order to deal with image effects, a phenomeno- 
logical fit is used. Given a good functional form /2(ri,r2) for P2 (invariant under simultaneous 
rotations of ri and Y2 as well as exchange of ri and r2), the nearest images are added in: 

/2(ri,r2)= £ /2(ri+iiiL,r2 + n2L). 

n\,n{e{-\fl,\} 

This function and its lattice integral /i(r) = ^''^'^) ^^"^ be simultaneously fit to p2 and pi 

using a nonlinear weighted least squares method. This allows the images to be subtracted off, 
giving p^ ~ p2 -/2 + /2 and p[ ~ pi -/i +/;, where /{(r) = / J3r2/^(r,r2). 

The so-called A ansatz for the static potential for interacting quarks [11, 12] was used as 
motivation for the functional form of f'j. We ultimately found that the following eleven parameter 
functional form gave a reasonably good fit: 

f^{r,,r2)=AgiruB,auO)g{r2,B,auO)g{r,C,a2,b2)e'"'^'^^ ^ ^ 



exp(— Ar) r> a 
ci — br — C2r^ r <a 



with g{r,A,a,b) 

where c 1^2 are given by the requirement that g and ^ are continuous at r = a. 



5. Lattice Calculations 



We used a mixed action scheme [13] with domain wall valence quarks on an asqtad sea, with 
mj, = 293(1) MeV and a = 0.1241(25) fm. This ensemble had 453 HYP smeared MILC gauge 
configurations [14], which have Nf = 2 + \ and volume 20^ x 64. Propagators were computed 
every 8 lattice units in the time direction, allowing for 8 measurements per gauge configuration, 
with source and sink separated by 8 lattice units. Wilson lines were computed using HYP smeared 
gauge links, and measurements were averaged over seven positions for the static quark: x = and 
the six nearest neighbors. 

For comparison, we also used a heavy quark mass, with ~ 940 MeV. Since the effects 
of dynamical sea quarks are negligible at that mass, we performed a calculation with K" = 0.153 
Wilson fermions on 200 configurations from the OSU_Q60a ensemble [15], which are 16^ x 32 
with quenched j3 = 6.00 Wilson action. From the static quark potential, this has a/ro = 0.186 [16]. 
Using ro = 0.47 fm, the lattice spacing is a = 0.088 fm. We used a source-sink separation of 1 1 
lattice units and averaged measurements over the two central timeslices. 

The functions /i 2 were fit to a restricted set of the lattice measurements pi .2- Three conditions 
were imposed to reduce the influence of the points most affected by images: r <^a for p\,r\ + r\< 
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Figure 2: Pi(r) without (left) and with (right) image corrections for the good diquark on the quenched 
m;i = 940 MeV ensemble 




Figure 3: P2(ri,r2) without (left) and with (right) image corrections for the good diquark on the quenched 
niji — 940 MeV ensemble, as a function of r with r _L R and = 0,4,6. 

lOOa^ for p2, and in both cases rimage > Ha, where nmage is the distance to the nearest periodic 
image of the static quark. Fits had psr degree of freedom ranging from 0.25 to 1.85. 

In the quenched good diquark case. Figure 2 shows the effect of image coiTcctions for pi, and 
here this procedure is quite successful, even extrapolating beyond the range included in the fit. 

For p2, the fit isn't as good as for pi , but it still works well. Figure 3 shows p2 with and without 
image corrections. The figure on the right looks cleaner for two reasons. First, the fit function is 
determined using a global fit, which allows for small deviation from the specified R and 6 to be 
compensated for. Second, image corrections have been applied, which are substantial for points 
distant from the origin. The end result is that the difference between the plotted point and the fit 
curve is equal to the difference between the raw data point and the fit for that point. 

6. Results and discussion 

As a check of how well the correlation function isolates the diquark from the effect of the 
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r (fm) r (fm) 

Figure 4: C{r\ , r2), as a function of r (in fm) with R = 0.2 fm (left) and R = 0.4 fm (right), with r _L R for 
the good and bad diquarks and the two pion masses. 

static quark, we compared different directions of r. Even at /? = 0.2 fm, C was independent of the 
direction of r, indicating that this correlation function works quite well. 

Finally, we can compare the systems. Fig. 4 shows the profile of the correlation function C 
at two fixed distances R from the static quark to the center of the diquark, and Fig. 5 shows the 
full dependence of C on r, at fixed R = 0.4 fm. The good diquark has a large positive correlation 
at small r that becomes negative at large r. The bad diquark has similar behavior with smaller 
magnitude. The difference between the good and bad diquarks is larger for the lighter pion mass, 
as expected from the quark mass dependence of the spin coupUng that splits good and bad diquarks. 
As R increases, both the correlation and the size of the positive region grow, although it is possible 
that some of this growth of C as /? increases may arise from the normalization of the correlation 
function. 

Our main conclusions are seen clearly in Fig. 5. The diquark correlations are highly indepen- 
dent of 6, indicating negligible polarization by the heavy quark, are much stronger in the good 
rather than the bad channel, and increase strongly with decreasing quark mass. Finally, it is impor- 
tant to note that the diquark radius is approximately 0.3 fm and the hadron half-density radius is 
also roughly 0.3 fm, so the diquark size is comparable to the hadron size. This is reminiscent of 
the size of Cooper pairs in nuclei, and argues against hadron models requiring point-like diquai^ks. 
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Figure 5: Continuous C(ri,r2) derived from the fit, as a function of r (in fm) with R = 0.4 fm. The two 
axes r|| and rx indicate directions of r parallel to and orthogonal to R, respectively. The color of the surface 
is discontinuous at C = 0. 
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